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Abstract 

We compute the force on a small neutral polarizable object moving at velocity v 
relative to a photon gas equilibrated at a temperature T We find a drag force linear 
in V. Its physical basis is identical to that in recent formulations of the dissipative 
component of the Casimir force. We estimate the strength of this universal Casimir 
drag force for different dielectric response functions and comment on its relevance 
in various contexts. 

PACS numbers: 72.20.-i, 72.20 Jv, 72.40.+W, 73.50.Pz 

The residual drag force on an AFM tip close to, but not in direct contact with a sub- 
strate in vacuo raises an important and fundamental question on the origin of non-contact 
friction PP. Other experimental techniques also are sensitive to non-contact friction j2]. 
Since it became clear that the Casimir effect, being par excellence a non-contact phe- 
nomenon, can lead to a dissipative drag (see ^ and references therein), it has become a 
primary focus of theoretical research Such Casimir dissipative drag occurs when elec- 
tromagnetic field fluctuations equilibrate in a specific reference frame, relative to which 
another system (e.g. a dielectric or a conducting body) is in uniform motion El E]- 
The difference between the frames of reference in the dissipative Casimir effect can be 
due to relative motion of different bodies, as for two conducting plates with relative mo- 
tion in the parallel direction, or for a neutral body moving relative to a conducting plate 
[3 El 13 Uni HI mi- In these cases the radiation equilibrates in one of the plates, and the 
friction depends upon the proximity of the other one. 

We show here that such friction can also come about when a single body moves 
relative to a thermal bath of the electromagnetic field excitations, such as those between 
the walls of an oven or in the cosmic microwave background. The friction has no position 
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dependence, i.e. it is spatially homogeneous. The consequence is a universal dissipative 
drag acting on all matter in relative motion with respect to a thermalized photon gas. To 
estimate the magnitude of this universal drag we evaluate it as a function of the dominant 
frequency of the electromagnetic response of the body for dielectrics and conductors. 
Consider the Lorentz force 

£f(^pE + jxBy (1) 

on a dielectric in a field E{r,t), B{r,t). The charge and the current densities in the 
dielectric are set by the polarization, P(r, t), such that p(r, t) = —V ■ P and j(r, t) = dtP, 
where P, j obviously obey the continuity equation dtp + V ■ j = 0. Assume that the 
response of the matter to the field is both linear as well as spatially local. Thus we write 

P(f, t) = eoJ Xe{t - t')E{f, t')dt\ (2) 

in MKS units, where the dielectric susceptibility Xe is dimensionless. In the frequency 
domain Xeijjj) is equal to (e(co') — l)/eo for condensed matter. For weakly interacting 
molecules or atoms Xe = PN<^m, with p^r the atomic or molecular number density of the 
medium and am the polarizability of the single molecule or atom. 

With the polarization proportional to the field, the force is bilinear in the field. This 
bilinearity holds for objects moving at arbitrary non-relativistic velocities relative to the 
frame of reference of the thermal bath. We require the thermal average of the force 
acting on a moving body in unbounded space filled with radiation at rest, working in the 
reference frame where the particle is instantaneously at rest and the photon gas moves at 
velocity v. The average force is obtained in terms of the thermal averages of the Fourier 
components of the field correlations. In unbounded space the Fourier components of the 
polarization and electric fields are 

P{k,iu),E{k,uj) = j dtj d'^fP{f,t),E{f,t) e-^(^-^^-^"*) (3) 

leading to the following form of the thermal average (..) of the Lorentz force 



(27r)4 J (27r) 
x{^{k' (^E{k, uj)E{k', cj')) - (k' ■ E{k, Lj)^ E{k\ J 



(4) 



where we took note of the fact that in empty space the electric field has no sources. A 
slight generalization of the standard expression (^21 Eq. 77.12) for the thermal average 
of the correlator of the electric field vectors yields 

(E,{k,u)E,ik\uj')^ = {2TiY5{k + k')5{uj + uj') {EiEj)^^^. (5) 
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with 



< EiEj >^^^-- 
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is the Bose occupation number for a photon distribution at temperature T, moving with 
velocity and (3 = (ksT)'^. Eq. IHlis obtained by generahzing the fluctuation-dissipation 
theorem to a translationally invariant system, and taking a Gibbs distribution correspond- 
ing to a photon gas moving with velocity v, just as done for excitations in superfluidity 
|13j . From Eq. IHlwe end up with the following average of the Lorentz force 
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Figure 1: Dependence of the logarithm of the relaxation time in years on the logarithm 
of the primary relaxation frequency (c^o in Eq- IT^ for three different temperatures, 
T = 300 K, T = 1000 K, T = 3000 K (upper, middle and lower curve). The frequency of 
the minimum is given by hwQ = 5.9694 /ebT. 



(7) 

Substitution of (4) in (5) with v = gives an integral over k that, by symmetry, is 
zero. To obtain a non-zero result it is convenient to first eliminate the delta functions by 
integrating over u and then expanding in powers of v (as done to obtain the normal fluid 
density in the theory of superfluidity). Integration over u then yields 



<F> = -i27rhcV J j^-^ kk Xe{ck) [l + 2n{ck,k)j + 

/d^k ^ / ^ \ 

-j—^kkxe{-ck)\l + 2n{-ck,k)j, (8) 
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where we have taken the dielectric to be homogeneous so the integral over its volume 
simply gives V. 

The same equation also can be obtained from energy loss considerations, following 
the arguments of Volokitin and Persson JU]. One starts from the dissipation rate of the 
energy in the rest frame of the fluctuating field, 



dW 



d^rij + pv) ■ E 



dWp 
dt 



-F-v. 



(9) 



Here we take into account that in the rest frame of the thermalized photon gas the total 
electric current is given via the linearized Lorentz form j — > j + pv. Apart from the heat 
production in the frame of the body, Eq. |Hlis reproduced immediately. 
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Figure 2: Dependence of the relaxation time in years on temperature. For the dielectric 
response (upper curve) we plot the minimal relaxation time obtained for the relaxation 
frequency that is proportional to the first Matsubara frequency. We take pm to be the 
density of water and the dimensionless measure of atomic polarizability ao = {N/V)eam ~ 
1. Numerically this value corresponds to a single particle of molecular polarizability 
coQ^m ~ 1-0 X 10~^° m'^ and a mass equal to the proton mass 1.67 x 10~^^ kg. For the 
metallic response (lower curve) we take the characteristic time Tq = ^ ~ 10^^^ s, well 
within the range of conductivities of simple metals. 

We now write, in a standard way, Xe{x) = Xe(^) +^Xe(2^)) noting that the real part is 
an even and the imaginary part is an odd function of the argument. The expression for 
the force can now be obtained in the following form 



<F>=-^^ I d^kkk 



x'e{ck) (^n{ck, k) — n{—ck, fc) j + 

+ i Xeick) (^l+n{ck,k) +n{-ck,k)^ . (10) 



Expanding the Bose function to the lowest order in velocity, we have n{ck, k)—n{—ck, k) = 
coth l/3hck + C(f2) and 1 + n{ck, k) + n{-ck, k) = Icsch^ l(3hck (^(3h{k ■ v)^ + 0{v^). 
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Placed in Eq. |H1 the A;-space integral over the term in x'ei^k) is zero, by symmetry. The 
remaining term then yields 

Performing the angular integral over d^k gives a factor of Att/3. Reverting again to uj by 
substituting k = u/c gives 

,p,^v,(m r^^i^. (12) 



This is the fundamental result of our paper. It states that the EM field fluctuations 
exert a drag, proportional in the lowest order to the velocity, on a particle that moves 
with respect to the frame of reference in which the EM field fluctuations are thermalized. 
Setting this force to Mw/r, where M is the total mass of the object and 1/r is the drag 
time, and using pm = M/V , yields the result that 

Regard this fundamental result in three different contexts: molecules, dielectric and con- 
ducting condensed matter. First assume that the dielectric response of the medium can 
be characterized by a single sharp absorption line at ujq. Because - is proportional to 
Xe(i^), obviously each of the absorption lines for a molecule or a dielectric will contribute 
additively to the integral in Eq. [T21 We set Xei}^) = Cio^{^/^o~^), where is a constant 
describing the strength of the response. This assumption gives 



/37mc^\ sinh^r^) / SnpMC^h^ \ sinh^ 



(14) 



after introducing x = \l3hjJo. In this form the relaxation time depends strongly on the 
absorption frequency ujq. It has a minimum at a temperature dependent frequency (see 
Fig. IH) that coincides with the minimum of the function f{x) = , at Xm = 2.98, 

where f{xm) = 0.137. Because f{x) is the square of the function it has a broad 

minimum with a quartic, rather than a quadratic, dependence upon the deviation from 
Xm- Taking this minimum at Xm into account, the smallest possible relaxation time can 
thus be obtained from the above equation in the form 

r = ro(|)', (15) 

where tqTq = ^ 37r/(a:„o j pm^^^ ^]-^jg minimum the absorption frequency hwo = 2xm ksT 

is proportional to the first Matsubara frequency. The temperature dependence of this 
minimal possible relaxation time is given in Fig. |21 
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A different formula for r is obtained for metals, which have constant conductivity at 



frequencies below the collision time of their charge- carriers. In this case Xei^) ~ 
Inserting this into Eq. ^1 the inverse relaxation time for drag now takes the form 



^ ' ^ duj o^, , (16) 



r \3t^Pmc^(^o J Jo sinh^(|/3/ia;) 
so that 

- = ro(|)', (17) 

with tq = ^ and Tq = ^^:^^r^- For most common metals the value of tq is between 
Xg-19 _ 10"!'' s. Taking its geometrical average we obtain the temperature dependence of 
the relaxation time as given in Fig. |21 

The times given in the figures are relatively long, corresponding to the general weakness 
of the Casimir interaction and the even lower strength of the dissipative Casimir effect. 
Two circumstances under which such long times might be observable are ovens and the 
cosmos. 

For ovens, on trapping molecules with unusually large polarizabilities it might be 
possible to observe a resonance in an ion or atom trap with a quality factor Q that 
is determined by the dissipative Casimir interaction. Small metal particles in a high 
temperature oven should also be susceptible to the effects of the Casimir drag. 

For the cosmos, it is believed that hydrogen atoms condensed from protons and elec- 
trons at about 3000 K, and that the coupling of radiation and matter due to Compton 
scattering becomes ineffective below this condensation temperature ^3]. However, as is 
clear from Fig. 2, it should not be difficult for molecules to remain coupled to the cosmic 
microwave background when the temperature was between about 300 K (and perhaps a 
bit less) and the 3000 K condensation temperature. This coupling could have an influence 
on the structure and anisotropies observed in recent experiments on the cosmic microwave 
background. It could also influence the behavior of molecules formed from the residue of 
novas and supernovas, and then subject to drag from a still-hot cosmic microwave (i.e., 
electromagnetic) background. 

The universal thermal drag calculated in this contribution is yet another important 
facet of the Casimir effect, which in this context appears to play a role not only in the 
static interactions between dielectrically inhomogeneous bodies but also in the dissipative 
dynamics of particles in homogeneous space. 
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